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The long-time, asymptotic state of rotating homogeneous shallow-water equations 
is investigated. Our analysis is based on long-time averaged rotating shallow-water 
equations describing interactions of large-scale, horizontal, two-dimensional motions 
with surface inertial-gravity waves field for a shallow, uniformly rotating fluid layer. 
These equations are obtained in two steps: first by introducing a Poincare/Kelvin 
linear propagator directly into classical shallow- water equations, then by averaging. 
The averaged equations describe interaction of wave fields with large-scale motions 
on time scales long compared to the time scale 1/ fo introduced by rotation (/o/2- 
angular velocity of background rotation). The present analysis is similar to the one 
presented by Waleffe (1991) for 3D Euler equations in a rotating frame. However, 
since three- wave interactions in rotating shallow-water equations are forbidden, the 
final equations describing the asymptotic state are simplified considerably. Special 
emphasis is given to a new conservation law found in the asymptotic state and 
de-coupling of the dynamics of the divergence free part of the velocity field. The 
possible rising of a de-coupled dynamics in the asymptotic state is also investigated 
for homogeneous turbulence subjected to a background rotation. In our analysis we 
use long-time expansion, where the velocity field is decomposed into the ‘slow man- 
ifold’ part (the manifold which is unaffected by the linear ‘rapid’ effects of rotation 
or the inertial waves) and a formal 3D disturbance. We derive the physical space 
version of the long-time averaged equations and consider an invariant, basis-free 
derivation. This formulation can be used to generalize Waleffe’s (1991) helical de- 
composition to viscous inhomogeneous flows (e.g. problems in cylindrical geometry 
with no-slip boundary conditions on the cylinder surface and homogeneous in the 
vertical direction). 


1. Introduction 

There are many important engineering (e.g, turbomachinery, ship propellers) 
and geophysical (geophysical turbulence) problems in which rotation significantly 
modify the turbulence properties of fluid flows. The interest in the effects of rota- 
tion is reflected in the large body of theoretical, experimental, and numerical work 
documenting them (Hopfinger 1989). As the effects axe both multi-fold and subtle, 
the development of models which account for the effects of rotation require an un- 
derstanding of the processes occurring in these flows. The objective of this work 
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is to examine the long-time evolution of solutions of shallow-water and 3D Euler 
equations in a rotating frame. We are especially interested in possible de-coupling 
of dynamics on slow manifold and in developing new conservation laws valid in the 
asymptotic state. 

The paper is organized as follows. The long-time averaged rotating shallow- water 
equations are presented in Section 2. Special emphasis is given to a new conserva- 
tion law found in the asymptotic state and to the de-coupling of dynamics of the 
divergence free part of the velocity field. Section 3 is devoted to homogeneous tur- 
bulence subjected to uniform rotation. In our analysis we use long-time expansion, 
where the velocity field is decomposed into the 'slow manifold’ part (the manifold 
which is unaffected by the linear ‘rapid’ effects of rotation or the inertial waves) and 
a formal 3D disturbance. We derive the physical space version of the long-time av- 
eraged equations and consider an invariant, basis-free derivation. This can be useful 
in generalizations of Waleffe (1991) helical decomposition to viscous inhomogeneous 
flows (e.g. problems in cylindrical geometry with no-slip boundary conditions on 
the cylinder surface and homogeneous in the vertical direction). Finally, the possi- 
ble rising of a de-coupled dynamics of the slow manifold is investigated. This study 
was motivated by the previous studies of Squires et al. (1993) and Cambon et al 
(1994). 

2. Rotating shallow- water equations 

It has been long realized that geophysical flows admit motions varying on differ- 
ent time scales. In many situations, the low frequency class contains the majority of 
energy and is therefore the class of main interest. However, long-time computation 
of this class of motions using unmodified Eulerian equations is prohibitive due to 
severe accuracy and time step restrictions (Browning et al 1990). The disparity 
of time scales leads to problems in the numerical solution of the equations because 
the Courant number is determined by the fastest time scale and therefore limits 
the time step which makes explicit solution impractical. In our approach collective 
contribution to the dynamics made by fast wave motions is accounted for by an 
averaging procedure. Elimination of “fast” time scales through the averaging pro- 
cedure leads to equations that more accurately and efficiently describe dynamics on 
long time scales. These equations are obtained in two steps: first by introducing 
a Poincare/Kelvin linear propagator directly into classical shallow- water equations, 
then by averaging. Averaged equations describe interaction of wave fields with 
large-scale motions on time scales long compared to a time scale 1/fo introduced by 
rotation (/ 0 /2-angular velocity of background rotation). These are fully nonlinear 
equations containing classical quasigeostrophic equation as a completely decoupled 
subsystem. Two other equations describing departure from quasigeostrophy are 
coupled to the latter equation. Special emphasis is given to a new conservation law 
found in the asymptotic state and de-coupling of dynamics of the divergence free 
part of the velocity field. 

The equations for shallow- water are: 

d 
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= -V • [VH] (2) 

where V and U are two-dimensional and in the (xi,x 2 )-plane, U is the velocity, 
g is the reduced acceleration of gravity (a constant), H is & function of Xj, x 2 , 
and t and is the free-surface height of the shallow layer of fluid, H 0 is a constant 
and is the mean depth of the layer, e 3 is the unit vector in x 3 , and fo is the 
Coriolis parameter. Equations (1) and (2) with boundary conditions describe the 
full system. The system (l)-(2) could also be thought of as describing a two- 
dimensional compressible gas, with H being the density. Note that ( 1 ) and (2) 
imply 

[^ + U • V](? = 0 (3) 

where Q is the potential vorticity 


Q = 


fo + fi 
H 


( 4 ) 


where 0 is the vorticity. Note that Eqs. (1), (3), and (4) are also a complete 
description of shallow- water motion. The fast time scale that we average over is 
l// 0 . Here we treat fo as constant, though later we will modify it to be a function 
of #2 or m polar coordinates a function of r. 

2. 1 Linear propagator 
We have for the linear propagator 

^ = - 3 Vif + / 0 Uxe 3 , H = -ff 0 VU. (5) 

Introducing D = V * U and SI = curl U • £ 3 , we obtain 

f + %—AD. f —A* <«> 

In the case of homogeneous flows we use Fourier representation (for simplicity we 
assume 2 tt periodicity in xj and x 2 ) 



k 


where k = (k u k 2 ), x = (xi,x 2 ). Then in Fourier space Eq. ( 6 ) becomes 


dDk i; |2rr , f o r t-v dH k 

= Hk + fo^lk, - 57 - = -foDk, 


Qt ~ dt J 

We supplement (7) with initial conditions 

Dk(0) = D°k, ft*(0) = nj, H k (0)=H° k . 


= -H 0 D k . 


(?) 


( 8 ) 



376 


A. Mahalov 


Solutions of the linear propagator Eqs. (7)-(8) are well-known (e.g. Pedlosky 1979); 
those corresponding to the free surface gravity waves represent the fast time scale. 
Defining = /o \/ 1 + | A: p where L 2 ^ = gHo j is the Rossby deformation 

radius, we have 


[W)\ /DJ\ 

\W)J \H°) 

( D ° k \ 

= exp(Lf) fij . 

\mJ 

Here A (£), A c (fc) and A , ( k ) are 3x3 matrices defined as follows 


(9) 


a (*o = 4 

<t>k 

a .(*) = 4 

<Pk 


( 0 0 0 \ 

0 gH 0 \k\ 2 —fog\k\ 2 , 

\0 -Hofo n ) 


( 0 /* g\k\ 2 \ 

\ ~fo 0 0 

\-Ho 0 0 / 


1 (*l 0 0 \ 

Ae(*) = -ry 0 f 2 f o9 \k\ 2 , 

\ 0 jy 0 /o ffi?o|^| 2 / 


( 10 ) 


£••2 The asymptotic limit of rotating homogeneous shallow-water equations 
We write classical rotating shallow-water equations (l)-(2) using the variables 
■D = V • U, Cl = cur/U • e 3 and H . Letting V <r = (D,Sl, H) (tr- transpose), these 
equations can be written symbolically in the form 


dV 

dt 


LV + B(V, V). 


( 11 ) 


Here L is the linear propagator operator corresponding to Poincare/Kelvin waves. 
We introduce this linear propagator directly into nonlinearity using the change 
of variables V = exp(Lt)v where v <r = (d,u>,h) and exp(Lt) is defined by (9). 
Equation (11) written in v variables has the form 

&v 

~Qt ~ B(<, v, v) (12) 

where 

B(£,v,v) = exp(—Lt)B(exp(Lt)v,exp(Lt)v). 

Equation (12) is explicitly time- dependent with rapidly varying coefficients. Using 
the standard averaging methods (e.g. Hale 1980), we average over the fast time 
scale of the gravity waves and obtain 


dv - . 

W = B (v, y ), 



v)ds. 


(13) 
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These equations are exact in the limit /i = / 0 T > 1 where /i is a non-dimensional 
parameter. We present explicit final form of the regularized (long-time averaged) 
Eqs. (13). The calculation leading to these equations is lengthy but straightfor- 
ward. Here we present only the final result for periodic case. We use the fol- 
lowing notations: m and k are wavevectors so that m = (mi,m 2 ), k = (ki,k 2 ); 
(mke 3 ) = (mxk)-e 3 = mik 2 - kim 2 ] k-x = kix 3 + k 2 x 2 and m k = mill +m 2 k 2 ; 

= f$+gH 0 \k\ 2 . The regularized (long-time averaged) equations are conveniently 
written in new variables 


A = d, B = fouj — gAh, C = — foh 


( 14 ) 


where A is the Laplace operator. Physical interpretation of these variables is as 
follows. C/Ho is quasigeostrophic potential vorticity, gAh/ fo is quasigeostrophic 
component of vorticity and, therefore, Bj fo is ageostrophic component of vorticity. 
In the quasigeostrophic case we have A = B = 0. Let 


A = £y n M„, b-jv—b,, c = ]ry n *c n . 


(15) 


We note that three-wave interactions are forbidden since the equation dzcpk ± Om — 
<f>„ (n = k + rrt ) has no solutions. Then the final equations describing the asymptotic 
state are much simpler than in the case of 3D Euler equations considered by Waleffe 
(1991). We obtain after averaging (Eq. (13) in A, B, C variables) 


dAn 

dt 

dB„ 


= C n - k (a(k,n - k)A k + P(k,n - k)B k ), 




4>l 


Here 


= V Cn-k{-<t>lP(k, n - k)A k + -rf ■ n - k)B k ), 
d ‘ * k 

it:no restriction 

/f x g(mke 3 )(2<t>l(m-k) + f£(\k\ 2 -\m\ 2 )) 

a{k ' m) = WMtZ ’ 

an , /o3(|m| 4 -|fc| 4 -||m| 2 |fc| 2 ) g(mke 3 ) 

«‘' m) mm. ■ y( ’ m) - 


(16) 


(17) 


are geometric parameters. Equations (16) are exact in the limit /i - /oT > 1. They 
have property that the third equation is completely decoupled from the first two 
equations. After the third equation is solved, the first two equations become linear 
equations in A and B. From the third equation in (16), it can be easily shown that 
C(t,x) obeys in physical space 

A c = 0 

Dt 


(18) 
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where is the advective derivative, based on the velocity v = e 3 x where 

V» = (1 — L 2 r A)~ 1 C/H 0 . Note that in the quasigeostrophic limit, ■tp is the stream 
function associated with the quasigeostrophic potential vorticity C / Hq, and v is the 
quasigeostrophic velocity. The above Eq. (18) implies that C and all of its moments 
are conserved because v is divergence free (a well-known fact within the context 
of two-dimensional incompressible turbulence). For a given geostrophic component 
C(t,x) found from (18), equations for ageostrophic components A and B have the 
form 
dA 

”= ^2 Cn-k{a(k,n - k)A k + (3(k,n- k)B k ), 

k:]k\=\n\ 

, 2 (19) 

= 5Z Cn-k{-4>l0(k, n - k)Ak -(- %«(*, n - k)B k ). 

*=I*I=W Vk 


dt 

dBn 

dt 


It can be easily shown that these equations have a new conservation law 


( 20 ) 


The conservation of M is proved by direct differentiation of (20), substitution of 
expressions for dA n j dt and dB n j dt from (19) and using properties of the geometric 
coefficients a(&,m) and j3(k,m) given by (17). Another approach for determining 
this conservation law would be to derive the long-time averaged equations from a 
variational principle. Existence of conserved quantity M allows us to control the 
size of ageostrophic component in terms of its initial value. One could use this 
conservation law to examine the Lyapunov stability conditions for the equilibrium 
solutions of the long-time averaged equations. We note that M is not conserved by 
the full rotating shallow-water Eqs. (l)-(2). 


On the asymptotic limit of 3D Euler equations in a rotating frame 

In a frame of reference rotating with constant angular velocity / 0 /2 about the x 3 
axis, the inviscid Euler equations have the form 

dU 

+ U • VU + /o JU = — Vp, V • U = 0. (21) 

Here U = (U, V , W) is the velocity field, p is modified pressure and J is the rotation 
matrix (JU = e 3 x U). It is convenient to write (21) in the form 

dU 

+ /o JU = U x cur/U - V(p + |U| 2 /2), V • U = 0. (22) 

We present a coordinate-free derivation of long-time averaged Euler equations de- 
scribing dynamics on time scales long compared to the time scale 1 / f 0 introduced 
by rotation. This can be useful in generalizations of Waleffe (1991) helical decompo- 
sition to viscous inhomogeneous flows (e.g. problems in cylindrical geometry with 
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no-slip boundary conditions on the cylinder surface and homogeneous in the vertical 
direction). 

We define the Leray projection P (projection on divergence free vector fields). 
Then Eq. (22) can be written in operator form as follows 

+ / 0 (PJ)U = P{U x curlV). (23) 

at 

We introduce a change of variables (canonical, preserving both energy and helicity) 

U = exp(-/ 0 PJ*)u. (24) 

This transformation is analogous to the van der Pol transformation widely used 
in the theory of nonlinear oscillators (e.g. Wiggins 1990). Equation (23) written in 
u- variables now becomes 

^ = exp(/oPJt)P{exp(-/ 0 PJOu x cur/(exp(-/ 0 PJ<)u)}. (25) 

at 

We note that the operator P J has the following representation (A = curl 2 is the 
Stokes operator) 

PJ = —A~ l -^—curl. (26) 

dx$ 

This representation in the periodic (homogeneous) case reduces to PJ(e'* r u*) = 
k x u* = e‘*' I j||iR(Jt)u* ; where u* is the Fourier coefficient corresponding 
to wavevector k = (fci,& 2 , ^ 3 ) R(fc)ut = k x ujt . The operator R(fc) becomes 
diagonal in the helical basis described in Waleffe (1991) (Eq. (3) in his paper). 
In the case of horizontally inhomogeneous flows a more general representation (26) 
should be used. From (26) we obtain 

exp(/ 0 PJt) = C(t) - A" ’ S (t)curl (27) 


and 

curl exp(/oPJt) = C (t)curl - A*S (t) (28) 


where 


B 9 

C (<) = cosh(/ 0 tA""&^-), S(t) = sinh(/otA~^^-). 


We note that C(t) is an even and S (t) is an odd operator-valued function. Us- 
ing (27), we have for the nonlinear term in (25) 


exp(/ 0 PJ<)P{U x curlV} = ( C(t ) - A~*S(t)cur/)P{U x curlV} = 

C(t){V x curlV} - A~*S(0[U,cur/U] - Vq (29) 
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since curl ( U x cur/U)= [U,cter/U]. Here [ , ] is a commutator of two vec- 

tor fields: [U, cur/U] = (U • V)cur/U — (cwr/XJ • V)U; q is a modified pressure. 
From (24), (27) and (28) we find that 

U = C(<)u + A~^ S(t)curlu, curlU = C(t)curlu + A^S(t)u. 

We substitute these expressions in (29). Then equation (25) becomes 

c?u 

-£■ = B(<, u, u) - Vg, V • u = 0. (30) 

Here B(t, u, u) is given by 

B (t, u, u) = C(<){U x curlU } - A“*S(t)[U, cur/U] = 
C(t){C(t)u x C(t)curlu) - A~iS(t)[C(t)u,C(t)curlu] + 
C(t){A~ ? S(t)curlu x A^S(f)u} - A" ^S(t)[A~^S(t)curIu, A*S(<)u]+(31) 
C(t){C(t)u x A*S(t)u} - A"*S(t)[C(<)u,A*S(<)u]+ 

C(t){A~?S(t)curlu x C(t)cur/u} - A-*S(t)[A-*S(t)cuWu,C(<)cur/u]. 

Separation of “fast” and “slow” oscillations in (31) is equivalent to a search for 
resonances which is a geometry dependent problem. One immediate simplification 
is that odd terms make no contribution to the time-averaged equations. Thus 
averaging of (30)-(31) reduces to averaging of four even terms 

^-'(O(C(f)u x C(<)ctxr/u}, C(t){A~iS(t)cur/u x A^S(t)u}, 
A-*S(t)[C(t)u,A*S(f)u], A-lS(t)[A-lS(t)curlu,C(t)curlu]. (32) 

It is convenient to work in the basis of the Stokes operator A. The operators C(<) 
and S(t) have a block-diagonal form in this basis. Since C(f) = cosh(/ 0 <A - ^ ^) 
and S(t) = sinh(/o<A - ^ g|j), they become simple oscillations in the basis of A. 
These oscillations have frequencies /o^- where \j are the eigenvalues of A and j 3 
is axial wavenumber. Then resonances in (32) are given by 

, h m 3 _ n 3 

vr k ^ 

In the C£ise of periodic (homogeneous) flows the helical basis formed by eigenfunc- 
tions of the curl operator acting in the space of divergence free vector fields is very 
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useful. The operator PJ ( jffrR(fc) in Fourier space) becomes diagonal in this basis. 
In this case we have Xj = \j\ where j = U 1 J 2 J 3 ) is the wavevector. Then averag- 
ing of (30) leads to equations studied by Waleffe (1991). In particular, he argued 
that in the asymptotic state of rotating homogeneous turbulence interactions are 
restricted to pairs of wavevectors satisfying resonance conditions (n = k + m) 


, *3 , m 3 _ n 3 

1*1 M H‘ 


(33) 


Now we consider homogeneous (periodic in three directions) flows and discuss 
equations describing the asymptotic limit of 3D Euler equations in a rotating frame. 
Both computations and experiments have noted an increase in integral length scales 
along the rotation axis relative to those in non-rotating turbulence. Increase in the 
integral length scales has been thought to be a prelude to a Taylor-Proudman re- 
organization to two-dimensional turbulence. It has been also suggested that this 
process is self-similar (Squires et al. 1993). The ‘slow manifold’ interactions (inter- 
actions which are unaffected by the linear ‘rapid’ effects of rotation or the inertial 
waves) are obtained in (33) with k 3 = m 3 = 0. These interactions are invariant 
under scalings of k 3 , m 3 in a trivial way. It is reasonable to conjecture that at 
zero Rossby number the only acting triads in the asymptotic state are those which 
are invariant under the scalings k 3 -> ffc 3 , m 3 -+ ■qm 3 where £ and 77 are arbi- 
trary real numbers. This is called the selection principle for triad interactions in 
the asymptotic state. It is inspired by the observations of Squires et al. (1993) on 
self-similar asymptotic states of rotating homogeneous turbulence. One can easily 
find solutions of the resonant Eq. (33) obeying the selection principle. These are 
{m 3 = 0, |n| = |Ar| } , {k 3 = 0, |n| = |m|} and {n 3 = 0, |fc| = |m|} with the signs 
chosen appropriately. Since these interactions are unaffected by scalings of the ver- 
tical wavenumber, they remain active even at zero Rossby number. It is possible to 
show that the interactions {n 3 = 0, \k\ = |m|} make no contribution to the averaged 
equations. Defining k± = (ki,k 2 ) to be the horizontal wavevector perpendicular to 
the rotation axis, interactions active in the asymptotic state can be written in the 
form {k 3 = m 3 = 0}, {m 3 — 0, |nx| = |*x|} and {*3 = 0, |nj_| = |rox|}. The form 
of these interaction manifolds suggests splitting the velocity field into two parts, 
each separately related to the horizontal or vertical velocities. We proceed to define 
splitting by first defining u> to be vertical (in the direction of rotation) component 
of the vorticity vector. Then we define splitting of the total velocity u as follows 

u = v + v where v = e 3 x VV* 


and tj> is related to u ; by the equation — — w. Applying averaging to the 

transformed Euler Eq. (30) (we note that only averaging of the even terms given 
by (32) is required) and using the selection principle, we obtain equations valid in 
the asymptotic state. Our analysis is similar to the one presented in Section 2.2 
where de-coupling of the divergence free part of the total shallow-water velocity 
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field was rigorously proved (see Eq. (18)). We find that in the asymptotic limit 
u>(#,x) obeys in physical space 

K" = 0 < 34 > 

where is the advective derivative, based on the velocity v = 63 x Vi/> where 

—Aip = u). This equation should be compared with Eq. (18). Thus we obtain 
de- coupling (in the asymptotic state) of the dynamics of v (two-component and 
horizontally divergence free). It satisfies classical 2D Euler equations (34). The 
de-coupling which is linked to the v part of the total velocity cam reflect the rise of 
an inverse energy cascade. Our future plans include confirmation of inverse cascade 
using the CTR database. Another avenue is to search for new conservation laws 
in the asymptotic limit of rotating homogeneous 3D Euler equations similar to the 
one we found for the long-time averaged shallow-water equations (Section 2.2, Eq. 
(20)). One possible approach for determining new conservation laws would be to 
derive the long-time averaged Euler equations from a variational principle. 
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